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Abstract 

In classical Kirchhoff flow, an ideal incompressible fluid flows past an obstacle and 
around a motionless wake bounded by free streamlines. Since 1869 it has been known 
that in principle, the two-dimensional Kirchhoff flow over a polygonal obstacle can be 
determined by constructing a conformal map onto a polygon in the log-hodograph plane. 
In practice, however, this idea has rarely been put to use except for very simple obstacles, 
because the conformal mapping problem has been too difficult. This paper presents a 
practical method for computing flows over arbitrary polygonal obstacles to high accuracy 
in a few seconds of computer time. We achieve this high speed and flexibility by working 
with a modified Schwarz-Christoffel integral that maps onto the flow region directly 
rather than onto the log-hodograph polygon. This integral and its associated parameter 
problem are treated numerically by methods developed earlier by Trefethen for standard 
Schwarz-Christoffel maps. 
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Notation 


z: physical space coordinate (region bounded by obstacle and free streamlines) 
w : velocity potential (slit plane) 
x = V2w/W +x. (upper half plane) 

£ = hodograph or conjugate-velodty (gearlike region or Riemann surface) 
fl= -log 5 (polygonal region or Riemann surface) 


1. Introduction 

Figure 1 shows the geometry we are concerned with. An ideal incompressible fluid 
in the complex z plane undergoes irrotational flow rightward past a solid obstacle I\ 



Figure 1. Geometry of the Kirchhoff flow problem. 

The complex velocity is denoted by v(z) and is normalized by v(») = l. The obstacle 
consists of n solid line segments r 4 = (z i _ 1 ,z i ), lsfcsn, bounded by vertices z k , Osfc<n, 
and -y^TT denotes the angle of r* counterclockwise from the horizontal. At an unspecified 
stagnation point z* along T, the flow divides between an upper part passing over z„ and a 
lower part passing under z 0 . At z 0 and z„, the fluid does not undergo infinite acceleration 
and turn 180° back on itself; rather, it continues smoothly towards z = +» by flowing 
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around a wake in which v^O. The curves of discontinuity separating the wake from the 
moving flow are two free streamlines with stream function equal to zero, which we label 
T_ and T+. The shapes of T_ and r + are unknown a priori, but are determined implicitly 
by the condition that all along both of them, |v(z)| must be constant and equal to 1. The 
physical origin of this condition is the requirement that the pressure should have the 
constant value p x throughout the wake, and be continuous across r_ and r+; |v| = l then 
follows from Bernoulli’s equation. 

Here then is our Kirchhoff flow problem: given an obstacle r, calculate the free 
streamlines r ± , stagnation point z„, and velocity field v(z) for a flow of the type described, 
together with associated numerical quantities such as lift and drag. 

Free-streamline flows of this kind have a long history. They were introduced by 
Helmholtz and Kirchhoff [8] in 1868-9 in an attempt to resolve D'Alembert * s paradox: in 
ideal potential flow (with no wake), the pressure forces around an object balance exactly 
and so the drag is zero [13]. Similar ideas apply also in the study of jets and cavities, 
where the free streamlines typically separate a liquid from a gas. We will say very little 
about the physical aspects of our problem, which is obviously idealized; for surveys of the 
large mathematical and physical literature of wakes, jets, and cavities, sec 
[2,5,7,9,10,13,19]. For some recent computational work in this area, see [1,3, 1H]. 

The Kirchhoff flow problem can be cast as a problem in complex analysis. Let £ 
denote the hodograph variable, which is simply the complex conjugate of velocity: 

ttO = v( z ). (l) 

Let G z denote the region of moving flow bounded by the solid boundary T and the 
unknown free streamlines r ± . Since the fluid is incompressible and irrotational, v (if 
interpreted as a real vector) is the gradient of a real velocity potential <j>(z) satisfying 
Laplace’s equation. Equivalently, { is the complex derivative of a complex velocity 
potential w(z) = 4>(z)+i\|»( z ), 

£ (0 - ( 2 ) 

where the stream function »|>(z) is conjugate to 4>. The function w(z) is analytic in G r , and 
maps G, conformally onto a slit plane G w = w(G t ), shown in Figure 2a, where the slit 
begins at the point w* = w(z,) at which the flow divides to go around the obstacle. 
Without loss of generality we can take iv,=0, so that G w is the slit plane C— R + . 
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Figure 2 . w and x domains. 

To simplify subsequent manip ulations, it is desirable to reduce G w to a half-plane. 
Let x be a new complex variable related to tv by 

tv = ^-(x—x.) 2 , x = V2tv/vy+x., (3) 

as indicated in Figure 2b, and write x k =x(w k ). Here W £ R and x • € (—1,1) are constants 
that will be chosen so that x 0 =-l and x„ = l. Now G x =x(G w ) is the upper half plane, 
with [-1,1] corresponding to T and with (-®,-l) and (1,») corresponding to T_ and T + , 
respectively. We write X k = x(r t ), X ± =x(r ± ). 

The classical hodograph method of solution for Kirchhoff flows begins by calculating 
a conformal map of G x onto the hodograph domain G { = £(G r ). What makes this possible 
is the fact that although G t is unkn own because of the free streamlines, Gj is (more or 
less) known. On the solid boundary, { has known argument, since the flow must be 
tangential. This argument is ir-y^ir or -* y 4 ir, depending on whether the point on the 
boundary lies downstream of z« in the direction of z 0 or z„, respectively: 

rir— forz^r*, x<x,, 

argi(z) = { _ 7i1T forz€r t , *>x.. ( 4a ) 

On the free streamlines, ( has known modulus: 

|{(z)| = l for z€T ± . (4b) 

Thus G { is a “gearlike” region bounded by circular arcs and subsets of rays passing 
through the origin. If we introduce the log-hodograph variable 


n(z) = -log5(z), 


( 5 ) 
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then the corresponding domain G fl is bounded by vertical and horizontal line segments. 
Therefore a conformal map of G x onto G n can be written explicitly in the form of a 
Schwarz-Christoffel integral. The classical solution method first calculates this map, 
thereby obtaining the analytic relationship between { and w, and then integrates (2) to get 
w and 5 as functions of z. 

In practice, the only flows that have ever been obtained in this way involve very 
-simple obstacles T, such as the flate plate considered by Kirchhoff and certain wedges 
[2,7,13]. As the number of sides increases the conformal mapping problem rapidly 
becomes too difficult for analytical solutions. The following are the principal reasons why 
this is true, and why even numerical solutions have not been carried very far. 

The first difficulty is that although the conformal map onto any polygon can be 
expressed by the Schwarz-Christoffel formula, this formula depends on accessory 
parameters or prevertices whose values must be determined numerically. This is the 
Schwarz-Christoffel parameter problem. In the past, researchers concerned with free- 
streamline flows have not had methods available for solving the parameter problem 
reliably. Our own interest in this project was motivated in part by the fact that the 
second author has recently developed a numerical method that does this, computing the 
map onto an arbitrary polygon typically in just a few seconds of computer time. This 
method has been implemented in a Fortran package called SCPACK [15,16]. 

A second difficulty is that even if the prevertices for the map fi(x) are known, the 
further integration of (2) must be performed to recover quantities of physical interest. As 
a result, Kirchhoff flow calculations based in the traditional way on the log-hodograph 
domain may be time-consuming, requiring on the order of minutes of computer time. 

The final and most serious difficulty is that except in the case of a very simple 
obstacle, G n is generally a polygonal Riemann surface rather than just a polygon, and 
contains slits and/or branch points of unknown dimensions and even of unknown 
topology. The occurrence of a Riemann surface is in itself not a serious problem, for the 
Schwarz-Christoffel formula can readily be modified to handle such domains [4]. But the 
presence of unknown dimensions and topology is more serious. To overcome this 
problem, one can formulate and solve a generalized Schwarz-Christoffel parameter problem 
[17] in which some of the conditions that determine £(z) involve dimensions in G : rather 
than Gq. We have successfully performed such calculations for certain geometries and 
hope to describe these results in a later paper. But it seems that this kind of calculation 
always requires careful attention to the details of the geometry of the hodograph domain. 
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As a result the construction of a single computer program to handle a wide range of 
obstacles by the hodograph method seems to be a difficult matter. 

In this paper we combine the numerical ideas of SCPACK with an analytical trick 
that appears in 54.1 of the book by Monakhov [10] (see also [18]) to calculate Kirchhoff 
flows over arbitrary polygonal obstacles much more efficiently than the classical 
hodograph method permits. The key idea is that since our flow problem has only a single 
pair of free streamlines meeting at infinity, z(w) can be written in the form of a modified 
Schwarz-Christoffel integral. By working with this integral directly, we dispense with all 
explicit consideration of hodograph domains, hence of Riemann surfaces and unknown 
slits and branch points, and we also avoid the need for a second integration. The result is 
a computer program that solves the problem of Figure 1 for an arbitrary obstacle to high 
accuracy typically in a matter of seconds, that is, at an expense comparable to that of 
finding the conformal map onto a closed polygon with the same number of vertices. 

Section 2 describes our modified Schwarz-Christoffel map for z(w), adapted from 
Monakhov. Section 3 outlines an efficient procedure for computing it numerically. 
Section 4 gives numerical results for some idealized problems. In later papers we will 
consider in greater detail certain particular configurations of physical interest, and we will 
also investigate models featuring a wake pressure p wake <p», which is essential if one 
wants quantitative agreement with laboratory data. 

Our Fortran package for Kirchho ff flows is a modification of SCPACK called 
KIR CHI. Machine-readable copies of SCPACK and KIRCH1 can be obtained by 
contacting the second author. 
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2. The modified Schwarz-Christoffel integral 

To put the modified Schwarz-Christoffel formula in context, we begin with a 
description of the standard Schwarz-Christoffel problem. Suppose P is a polygonal region 
in the complex z-plane with vertices z k , 1 sjts n , and sides T k = (z k - u z k ) oriented at angles 
y k Tt counterclockwise from the real axis, and define the external angle parameters 3* by 

Pi = *yi+i“7t. l^rss/J-1 (6) 

a™ 1 Pn = 7i~7„+2. (For convenience we write z 0 =z n , z„ +1 =z lf etc.) Let z=z(x) be a 
conformal map of the upper half complex x plane onto P, with z(x=»)gr 1 , and let 
*i=*( r i) be the interval (x^.x*) bounded by the prevertices x k .^x(z k . x ) and x k =x(z k ), 
where x, < x 2 < • • • <x„. Then arg^(x) is a known function on all of R which has 
constant value on each X k , and jumps by 3*1* at x k : 

= 7i^ for x € (x„ , °o), (7a) 

= 3 i 17 at x=x k , 1 stSfi, (7b) 

The basis of the Schwarz-Christoffel formula is the fact that it is easy to write down the 
function ^ determined by these conditions. Let g k be defined by 

&(*) = ( 8 ) 

with the branch chosen so that g k (x) is positive for x>x k . Then g k has constant argument 

on R except for a jump by 3*^ at x k . To be precise, it maps lmx>0 onto the wedge 
bounded by the rays e ,p,1T R + and R + , as shown in Figure 3a. It follows that £ can be 
written as a product of these wedge maps, 

£ « - a ( 9 ) 

i-i i-i 

for some A > 0. The Schwarz-Christoffel formula is the integral of this , 

z(x) = C +aAY n (*' -x k )~*'dx ' , (10) 

k— 1 

where C is a complex constant. 

Returning to the Kirchhoff flow problem, consider now the function z(w) that maps 
the slit plane G w of Figure 2a onto the flow domain G t of Figure 1. Instead of knowing 

for ^ x€R » ( 2 )*( 4 ) we now know arg^ for x€[— 1,1] and |^| elsewhere. 
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Figure 3. Mapping properties of the individual factors g k and h k in the 
Schwarz-Christoffel and modified Schwarz-Christoffel integrands, respectively. 


Defining p t again by (6) and setting P« = 1, we can write 



= In* 

for x=x„, 

(Ha) 


= 

a\x=x k , lsfcsn-1 and *=“*”, 

(lib) 

1 dz 1 
| dw \ 

= 1 

for x£X ± , i.e. |jc| > 1, 

(11c) 

dz 

*8*7 

= 0 

at x = «. 

(lid) 


Thus the Kirchhoff flow problem is a modification of the Schwarz-Christoffel problem (7) 
in which a constant-modulus condition rather than a constant-argument condition is 
applied over part of the boundary. As soon as one formulates the problem in these terms 
it becomes evident that here again, ■— can be written as a product. Let h k be defined by 



h k (x) = 


( 12 ) 
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and h, similarly, where the branch is taken so that h k (x) is positive for x € (x k , 1). 
Obviously h k has a singularity like that of g k at x 4> plus additional sin gnlaritiV* at x= ±1. 

In fact, h k maps lmx>0 onto the pie-slice-shaped region bounded by the ray the 

ray R + , and the circle |z| = l, as shown in Figure 3b. (If p 4 >0, the region becomes an 
inverted pie slice with sides meeting at «.) Because every h k has modulus 1 outside 
[-1,1], we can write ^ as 

t - f - < v nw (u) 

Here and below, JJ denotes the product over *=1,2, . . . ,n-l and *=“*”, and Y is 
* * 

the analogous sum. 

By construction, the function on the right in (13) satisfies all of the conditions (11) 
except possibly (lid). We must choose x. so that (lid) is satisfied too. To do this, note 
that by an easy computation 

argA^oo) = -p 4 cos _1 (— x*), 
with the analogous formula for A*, and therefore 

“8 (^IIW»)] = 7 n '*r-2Pi C0S " 1 (“*4)- 
Condition (lid) therefore amounts to 

= 0 , 

* 

that is, 



The Kirchhoff flow can now be written as the integral of (13), 

z(w) = C+ e^f (15) 

For this formula to be usable we want to integrate with respect to x rather than w . By 
(3), dw' can be replaced by W(x'-x»)dx', and we get 


z(x) = C + wf—^-dx' = C + (x' -x,) X\h k (x')dx' 


(16) 
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in direct analogy to (10). 

The factor ( x'-x .) in (16) cancels another factor (x'-x*) -1 hidden there in h,. 
Consequently z(x) has singularities at x k , 0 but not at x*. We can write the 
integral out in full as follows: 


z(x) = C + We‘ 1, ' r f (l-x.x+'V (1— x 2 )^— x, 2 ^ J| 


x—x k 

l-x^+Va-^a-jD 


-P* 

dx\ (17) 


This is essentially eq. (5) of p. 185 of [ 10 ], except that Monakhov unjustifiably assumes 
x. = 0 and hence w> 0 = w„. 
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3. Numerical solution procedure 

The last section reduced the Kirchhoff flow problem of Figure 1 to the modified 
Schwarz-Christoffel integral (16). As always in Schwarz- Christoff el mapping, before one 
can make use of this integral one must solve a parameter problem in order that not only 
the angles but also the lengths of sides in G t come out right. In this case the parameter 
problem involves n-1 unknown prevertices {xj: 

n-1 unknowns: x lf x 2 , . . . ,x„_i satisfying — l<x : < • • • <x„_ 1 <l. (18) 

The correct values for {xj will be determined by an iterative process. Suppose that at 
some step of this iteration a set of estimates {xj is available. Then we first calculate a 
corresponding value x* from (14), and a value W from the side-length condition 

* i z i~ z oi 

r * 

X 0 

derived from (16). This leaves n— 1 further side-length conditions to be satisfied: 

- T *+l 

n-1 equations: \Wf(x’— x»)\\h k (x')dx'\ = \z k+l —z,\, lsJlSn-1. (19) 

Thus the count is right, and under suitable additional hypotheses, the existence of a 
unique solution to (18)-(19) can presumably be proved. . For general {xj these equalities 
will not hold, and the errors in them will be used to devise a new guess {x fc } for the next 
iterate. 

We now sketch how we carry out this process numerically in the Fortran package 
KIR CHI. Most of the ideas are adapted from SCPACK and arc discussed more fully in 
[15]. 

Solution of nonlinear system of equations. There is little reason to write one’s own 
program to solve (18) -(19); excellent robust programs for this purpose already exist in the 
public domain. Several of these arc based on the “hybrid” methods developed by M. J. 
D. Powell, which combine a steepest descent algorithm in early stages with a quasi- 
Newton algorithm as the solution is approached. One well-known program of this kind is 
the HYBRD1 code in the MINPACK library from Argonne National Laboratory, U.S.A. 
All of our own work has used instead Powell’s code NS01A [11], which can also be found 
in the Harwell Subroutine Library. Although one could in principle compute the Jacobian 
matrix for (18)-(19) exactly, NS01A achieves superlinear convergence while requiring 
function values only. Beginning with the trivial initial guess of equally spaced prevertices, 
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KIRCH1 typically converges to machine accuracy in around 4 n iterations. 

Change of variables to eliminate constraints. The parameter problem as written is 
constrained by the prevertex ordering conditions (18). If these constraints are ignored, 
TTyaningW* iterates will be generated and the correct solution will generally not be found. 
However, it is a simple matter to eliminate the constraints by adapting the change of 
variables introduced in SCPACK and also proposed by Reppe [12]. Define 

y = log^l, l^k^n-1. (20) 

x k+1 -x k 

Then there is a one-to-one correspondence between sets of unconstrained parameters 
{yjeir -1 and sets of constrained parameters {xj satisfying (18). To work in these 
variables, one can simply write the function evaluation subroutine called by NS01A so 
that it takes as input {yj instead of {xj. The effect of this is that in our experience so 
far, KIRCH1 has converged to the correct solution for every obstacle T attempted. 

Compound Ganss-Jacobi Integration. Obviously (16) must be evaluated 
numerically except in the most trivial cases. Like the standard Schwarz-Christoffcl 
integrand of (10), the integrand in (16) has the form A(x)(x-x i )" p * for some analytic 
function h near each prevertex x k , lstjfc=sn-l, and the numerical integration procedure 
must take this singularity into account or it will be hopelessly inefficient. In KIRCH1 we 
apply the appropriate Gauss- Jacobi quadrature formula, whose nodes and weights are 
determined by calling the program GAUSSQ by Golub and Welsch [6], a version of which 
is also available in the NAG Subroutine library. However, in the development of 
SCPACK it was observed that Gauss-Jacobi quadrature alone is not enough to ensure 
accurate integrals, because the exponentially large crowding factors common in conformal 
mapping often lead to one prevertex x } being so near another one x k that the associated 
singularity strongly affects intervals of integration ending at x k (see also [12]). The 
solution devised there was to divide the interval into subintervals on which Gauss-Jacobi 
or pure Gauss (-Legendre) rules are applied, with the lengths of these subintervals chosen 
dynamically in such a way that none is ever longer than the distance to the nearest 
singularity [15]. This is what is meant by compound Gauss-Jacobi quadrature. We 
recommend the same procedure for Kirchhoff flow computations. It has proved highly 
effective in KIR CHI, where we consistently obtain integrals accurate to around d digits 
when 2d is taken as the standard number of quadrature points per interval. 



- 12 - 


Numerical integration near the separation points. At the separation points x 0 and 
x ni the integrand of (16) has a singularity that is not simply of Gauss-Jacobi type, for 
although the boundary of G t does not turn a comer at these points, it changes abruptly 
from a straight segment to a smooth curve. Again the singularity must be treated 
properly if the integration is to be efficient. To determine its form, note tha t at x 0 , say, 
the boundary of the log-hodograph domain G^ defined by (5) consists of an intersection of 
a vertical and a horizontal straight line segment. Therefore ft(x) has the form 
n(x) = fl(*o) + A(x)V x-x 0 near x 0 with h analytic. It follows that the integrand of (16), 
say H(x), has a singularity of the same type as at x 0 , i.e. H(x) = A(VS^g for some 

new analytic function h. To integrate this, introduce a new variable y = Vx'-x 0 . Then 
the integral becomes 

x Vr-1 0 

fH(x’)dx = f lyHtf+x^dy, (21) 

*o 0 

and the integrand on the right is analytic at y = 0, so it can be treated by Gauss-Legendre 
quadrature. We apply this rule in the usual compound way. 

Evaluation of the inverse map. Equation (16) gives z as a function of x, which is 
what is needed for solving the parameter problem or producing plots of streamlines and 
equipotential lines as in the next section. To determine the potential or velocity at a given 
point in space, however, one needs to evaluate the inverse map x(z). The obvious 
approach to this is to solve the equation z(x)=z iteratively for x by Newton’s method. 
Since — is known exactly as the integrand of (16), this iteration can be carried out very 

efficiently. In practice any reasonable initial guess for x typically leads to convergence in 
3 or 4 iterations. 

Computation of drag and lift. By Bernoulli’s equation, the pressure at a point z € T 
acts in the normal direction pointing into the fluid and is equal to to y |v(z)| 2 = if we 

assume the fluid has density 1 and take the pressure for v=0 as zero. Therefore the total 
force on a segment of r bounded by z a and z b is 

F *T m = ±JX(w)dw = +ff(x-x.);(x)dx = lKf( x -x,)X(x)dx. (22) 

Now as pointed out in sightly different formulations by Levi-Civita and by Scbi ffman (see 
[5], p. 370 and p. 350), this result has a remarkable interpretation. Since {(*) maps 
(-oo, -i) onto an arc of the unit circle, it can be analytically continued by reflection to a 
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function £ in the lower half x-plane defined by the formula £(x) = l/£(x) for Imx^O. 
Thus (22) is equivalent to 

f2F = ^2~S^jrr~ dx - ( 23 ) 

By (16), if 5 were £, this would be a formula for -^(z b -z a ). In other words to determine 

Fffi, we have reflected r across the free streamline T_ into a new obstacle t , and each 
force along T corresponds to a distance along f : 

FSF = ( 24 ) 


This interpretation of the forces on r suggests immediately how to compute them 
numerically: one calculates the dimensions of t by the usual compound Gauss-Jacobi 
quadrature procedure, makin g use of the branch £ instead of £ and of Gauss-Jacobi 
formulas based on exponents 0* = -0*. 

In the wake, the pressure has the constant value corresponding to |v x | = 1, so the 
forces are given by 


*3* - i'U = 




(25) 


where we have assumed that the wake lies to the left as T is traversed from z a to z b . 
From (23)-(25) it follows that the total force on the obstacle from both sides is 

F = F^+Fffi = j(£ n -z n ), (26) 


since z 0 = z 0 . It is customary to resolve this number into its components. 


F - ^diag+^iift ~ ( C d +C l )L , (27) 

where the drag and lift coefficients C D and C L are defined in terms of some reference 
dimension L of T, such as its total height or length. 

Together, £ and £ define a single-valued analytic function £ on C\[-l,l], and the 
force can equivalently be obtained by integrating (25) clockwise around any contour 
enclosing [-1,1] in this domain: 


_ tW r (*-*.) , 
F 2 * £(x) *- 


(28) 


Clearly this number is equal to 2ir»‘ times an appropriate residue at but in practice it is 
just as convenient to obtain F from (26) by computing integrals for and z„. 
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4. Computed results 

To test the accuracy of KIRCH1, we have compared lift and drag coefficients 
computed by it with various values that arc given in the literature. First we considered 
three geometries whose exact solutions are known: the flat plate due to Kirchhoff in 1869 
[5, p. 329], the inclined plate due to Rayleigh in 1876 [5, p. 330], and the symmetric 
wedge due to Bobyleff in 1881 [9, p. 104]. In these cases KIR CHI reproduces the exact 
values C D and C L to many digits. Then we considered certain geometries for which 
numerical solutions have been published: an inclined plate with separation from the back 
face studied by Chaplygin and Lavrentiev in 1933 and by Sekerzh-Zenkovich in 1934 [7], 
an asymmetrical wedge studied by C. C. Lin in 1960 [20], a symmetrical 4-piece wedge 
studied by Wu and Wang in 1964 [20], and a plate with spoiler studied previously by the 
first author [3]. (The values tablulated under the heading “LS” in [3] are in error, and 
were replaced for this comparison by Herat’s previously announced corrections.) We also 
considered a circular arc of half-angle 55°, studied by Brodetsky in 1923 and Schmieden in 
1929 [13], which we approximated by inscribed polygons. In all of these cases KIRCH1 
reproduces the published values C D and C L up to small discrepancies which we attribute to 
the published sources, except that the numbers of Sekerzh-Zenkovich reported in [7] 
appear to be wrong. 

On the basis of these tests and other evidence, we believe that KIRCH1 can reliably 
compute flows to arbitrary accuracy over arbitrary polygonal obstacles with up to one or 
two dozen vertices. The number of correct digits obtained increases roughly linearly with 
the number of Gauss-Jacobi quadrature points per interval, and therefore the accuracy can 
be doubled by roughly doubling the computation time. 

We will now present some new Kirchhoff flow calculations of our own, summarized 
in Figures 4-9. Our purpose in presenting these examples is first, to demonstrate visually 
that the ideas described above really work, and second, to record some further numbers 
for comparison with future experiments. Examples of more direct physical interest will be 
considered in later papers. 

Each figure shows the obstacle T together with a system of streamlines at invervals 
AtJ/ = .1 and equipotential lines at intervals A4> = .2. These curves are obtained by 
mapping a rectilinear grid in G w conformally into G,. The first thing to notice in looking 
at the figures is that evidently the Kirchhoff flow problem has indeed been solved - for 
along the free streamlines, the equipotential lines are in each case evenly spaced, 
signifying constant flow speed. Moreover, since the successive lines are separated by 
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distances of .2, the speed has the correct value 1. By contrast, note that in the examples 
with large drags especially, the separation between adjacent streamlines or equipotential 
lines in the inflow end of the plots is visibly greater. 

The additional shape shown inside each wake (in Figs. 6 and 8, below the wake) is 
the reflected boundary £ described in the last section. Long sides in £ correspond to high 
flow speeds. The segments of £ that come out especially small are those along which the 
fluid is nearly stagnant. 

Each figure lists the computed drag and lift coefficients and parameter W, all of 
them probably accurate to the six digits given. The number of steps in the iterative 
solution of the parameter problem by NS01A is also listed. The total amount of work for 
solving the parameter problem scales roughly as n 2 times this number of steps, and for 
drawing the plot, as n 2 times a large constant [15]. A convenient and machine- 
independent way to measure the computer times for these tasks is to count the total 
number of complex logarithms calculated in all products (17) during the computation, for 
these calculations turn out to dominate the total computer time. In each figure the first 
logarithm count listed corresponds to the solution of the parameter problem, and the 
second to the construction of the plot. The counts are approximate. 

It remains to state the dimensions and reference lengths L of the various obstacles, 
and to make a few comments on each. 

Figure 4: flat plate. This is the problem treated by Kirchhoff, whose solution can 
be found in several of our references, e.g. [5], The plate has length L = 1. The exact 
values for C D and W are 2 'tt/(4+it) and 2/(4+ ir), respectively. 

Figure 5: wedge. Here the plate still has length L = 1 but has been bent at the 
middle, with the lower half inclined at an angle 45° and the upper half at 30°. The 
computed stagnation point lies on the lower face at a distance .017348 from the vertex. 
Even simple wedge problems of this kind cannot be solved analytically. 

Figure 6: plate with separation from rear face. This obstacle, motivated by the 
model proposed by Chaplygin and Lavrentiev mentioned above [7], consists of a plate of 
length L = 1 inclined at angle 30° that bends back 180° at the leading edge into another 
plate of length 1/2. In other words, it is an inclined plate with separation prescribed at 
the middle of the back face. Note that the sharp edge in T maps to a broad channel 
attending to « in £ . The stagnation point lies at a distance .193308 from the leading 
edge. 
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Figure 7: plate with spoiler. Here a spoiler of length .2 at angle 45 ° has been added 
to the last obstacle, forming the geometry considered previously by the first author [ 3 ] 
and also by Bassanini [1]. The stagnation point moves up to a distance .090124 from the 
leading edge, and the singularity there weakens considerably. 

Figure 8: equilateral triangle. This obstacle is a symmetrical equilateral t riangle 
with side lengths L = 1/2, having four vertices all told since the two separation points are 
mathematically distinct. The flow crosses over itself, hence is nonphysical. We give this 
example to emphasize that nothing in our formlulation requires that the Kirchhnff flow, 
or even the obstacle itself, be embeddable in the plane. We have also computed some 
extremely nonphysical flows over various more exotic obstacles, but there is no space to 
present them here. 

Figure 9: charm bracelet. Finally, we include this 13-gon, whose aerodynamic 
importance is limited, to emphasize that the methods described in this paper work for 
arbitrary polygons. The three straight legs, inclined at angles 45°, 90°, and 30°, each have 
length 1 / 2 , and they meet each other at distances .15 and .1 above the square and the 
triangle, respectively, both of which have have side lengths .15. The reference length is 
L = 1. Note that the sections £^£3 and £ fi -£ 10 of £ are so small, since the fluid is nearly 
stagnant there, that they are scarcely visible in the plot. 
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